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Abstract 

We investigate the possibility of deriving analytical formulas for the 15-dimensional separable 
volumes, in terms of any of a number of metrics of interest (Hilbert-Schmidt [HS], Bures,...), of 
the two-qubit (four-level) systems. This would appear to require 15-fold symbolic integrations 
over a complicated convex body (defined by both separability and feasibility constraints). The 
associated 15-dimensional integrands — in terms of the Tilma-Byrd-Sudarshan Euler-angle-based 
parameterization of the 4x4 density matrices p (J. Phys. A 35 [2002], 10445) — would be 
the products of 12-dimensional Haar measure pHaar {common to each metric) and 3-dimensional 
measures ^metric {specific to each metric) over the Sd-simplex formed by the four eigenvalues of 
p. We attempt here to estimate/determine the 3-dimensional integrands (the products of the 
various [known] PmetricS and an unknown symmetric weighting function W) remaining after the 
(putative) 12-fold integration of pHaar over the twelve Euler angles. We do this by first fitting W , 
so that the conjectured HS separable volumes and hyperareas {Phys. Rev. A 71 [2005], 052319; cf. 
quant-ph/0609006 1 are reproduced. We further evaluate a number of possible such choices of W 



by seeing how well they also yield the conjectured separable volumes for the Bures, Kubo-Mori, 
Wigner-Yanase and (arithmetic) average monotone metrics and the conjectured separable Bures 
hyperarea {J. Geom. Phys. 53 [2005], 74, Table VI). We, in fact, find two such exact (rather 
similar) choices for W that give these five conjectured (non-HS) values all within 5%. In addition 
to the above-mentioned Euler angle parameterization of />, we make extensive use of the Bloore 
parameterization (J. Phys. A 9 [1976], 2059) in a companion set of two-qubit separability analyses. 

PACS numbers: Valid PACS 03.67.-a, 2.40.Dr, 02.40.Ft, 02.30.Uu 
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I. INTRODUCTION 

In a pair of major, skillful papers, making use of the theory of random matrices 
Sommers and Zyczkowski were able to derive explicit formulas for the volumes occupied 
by the d = {n^ — l)-dimensional convex set of n x n (complex) density matrices (as well 
as the d = ^""^If""*"^-* -dimensional convex set of real n x n density matrices), both in terms 

n 

of the Hilbert-Schmidt (HS) metric — inducing the fiat, Euclidean geometry — and 
the Bures metric ^ (cf. ^). (These results are also more lately discussed in the highly 
comprehensive new text of Bengtsson and Zyczkowski 0, chap. 14].) Of course, it would 
be of obvious considerable quantum-information-theoretic interest in the cases that n is a 
composite number, to also obtain HS and Bures volume formulas restricted to those states 

n 

that are separable — the sum of product states — in terms of some factorization of n |6|. 
Then, by taking ratios — employing these Sommers-Zyczkowski results — one would obtain 
corresponding separability probabilities. 

In particular, again for the 15-dimensional complex case, n = 4 = 2 x 2, numerical 
evidence has been adduced that the Bures volume of separable states is (quite elegantly) 



2-i5( VO) ~ 4.2136 ■ 10-6 fl, Table VI] and the HS volume (5\/3)-^ ^ 2.73707-10-^ js!, eq. 



(41)]. Then, taking ratios (using the Sommers-Zyczkowski results we have the derived 

conjectures that the Bures separability probability is ^^'^Q^^^^) ~ 0.0733389 and the HS one, 
considerably larger, ^^'^''^^4^6^^'^ ~ 0.242379 ^ eq. (43), but misprinted as 5^ not 5"^ there]. 
(Szarek, Bengtsson and Zyczkowski — motivated by the numerical findings of 0, 0] — have 
recently formally demonstrated "that the probability to find a random state to be separable 
equals 2 times the probability to find a random boundary state to be separable, provided 
the random states are generated uniformly with respect to the Hilbert-Schmidt (Euclidean) 
distance. An analogous property holds for the set of positive-partial-transpose states for 
an arbitrary bipartite system" [ij (cf. ^|). ("Since our reasoning hinges directly on the 
Euclidean geometry, it does not allow to predict any values of analogous ratios computed 
with respect to Bures measure, nor other measures" [l^ p. L125].) These three authors 



also noted 



lOl . p. L125] that "one could try to obtain similar results for a general class 



of multipartite systems". In this latter vein, numerical analyses of ours give some [but 
certainly not yet conclusive] indication that for the t/iree-qubit triseparable states, there is 
an analogous probability ratio of 6 — rather than 2.) 
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However, the analytical derivation of (coniecturally) exact formulas for these HS and 



Bures (as well as other, such as the Kubo-Mori 1^ and Wigner-Yanase 



separable 



volumes still appears to be quite remote (cf. 3]) — the only such progress to report so far 
being certain exact formulas when the number of dimensions of the 15-dimensional space 
of 4 X 4 density matrices has been severely curtailed (nullifying or holding constant most 

(cf. Q]). Most notably, in this research direction, 



of the 15 parameters) to d < 3 

n 

in [16|, Fig. 11], we were able to find a highly interesting/intricate (one-dimensional) con- 
tinuum (— oo < (3 < oo) of two-dimensional (the associated parameters being bi, the mean, 
and (Tq, the variance of the Bell-CHSH observable) HS separability probabilities, in which 
the golden ratio was featured, among other items. (The associated HS volume element 
— 32/3(i+/3) dpdbqdag — is independent of bi and in this three-dimensional scenario. Ex- 
tensions to higher-dimensional scenarios d > 3 appear problematical, though.) Further, in 
isl ]. building upon work of Jakobczyk and Siennicki |l^, we obtained a remarkably wide- 
ranging variety of exact HS separability {n = 4, 6) and PPT (positive partial transpose) 
(n = 8, 9, 10) probabilities based on two-dimensional sections of sets of (generalized) Bloch 
vectors corresponding to n x n density matrices. 

In this paper we are able to report some additional progress in these directions. We 
obtain exact formulas for certain = 4, n = 4 scenarios and upper bounds for d = 7 and 
d = 9 instances. (Nevertheless, the full d = 9 and/or d = 15, n = A real and complex 
scenarios still appear quite daunting — due to the numerous separability constraints at 
work, some being active [binding] in certain regions and in complementary regions, inactive 
[nonbinding] . "The geometry of the 15-dimensional set of separable states of two qubits is 
not easy to describe" p. L125].) 

To proceed initially fsecs. ITTllllllll V|) . we employ the (quite simple) form of parameteriza- 



tion of the density matrices put forth by Bloore j2 
are several other possible parametrization s |2lL l22 



sorne thirty years ago. (Of course, there 
a number of which 



some thirty years aj 

y,y,y,y,Q, 



we have also utilized in various studies j2Sl, |22| to estimate volumes of separable states. Our 
greatest progress at this stage, in terms of increasing dimensionality, has been achieved with 
the Bloore parameterization — due to a certain computationally attractive feature of it, 
allowing us to decouple diagonal and non-diagonal parameters — as detailed shortly below.) 

In our final (quite differently structured) series of analyses (sec. |Vj), though, we em- 
ploy not the Bloore parameterization, but the Euler-angle-based one of Tilma, Byrd and 
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Sudarshan 



23|. Our motivation here is to bypass/circumvent the necessity of the puta- 
tively achieveable, but computationally daunting first twelve steps (over the twelve Euler 
angles) of a 15-fold integration, ("we would like to derive a subset of the ranges of the 
Euler angle parameters... dividing the 15-parameter space into entangled and separable sub- 
sets. Unfortunately, due to the complicated nature of the parameterization, both numerical 
and symbolic calculations of the eigenvalues of the partial transpose... have become com- 



putationally intractable using standard mathematical software" j23|, p. 10453].) We seek 
to find the t/iree-dimensional (weighting) function iy(Ai, A2, A3, A4) of the /owr eigenvalues 
(A4 = 1 — Ai — A2 — A3) which would yield our various conjectured 15-dimensional separa- 
ble volumes and 14-dimensional separable hyperareas for the two-qubit systems, ^^y]- In 
sec. IVll we also apply this approach exploratorily to the t?, = 6, = 35 case of quhii-qutrit 
pairs, seeking to find a suitable 5-dimensional weighting function VF„=6(Ai, A2, A3, A4, A5, Ar). 

Let us also point out to the reader, our even more recent companion-type study j30|, in 
which solely for the (relatively simple) Hilbert-Schmidt (non-monotone) metric, are we able 
(making use of the Bloore parameterization) to reduce the two-qubit (both real and complex) 
separable volume determination problems to those of finding one-dimensional weighting 
functions (which turned out to be well-approximated by certain incomplete beta functions 
— functions of the ratio of diagonal entries of the 4x4 density matrices p). 

P22P33 ^ J r J 



II. BLOORE PARAMETERIZATION OF THE DENSITY MATRICES 

The main presentation of Bloore [2^ was made in terms of the 3 x 3 (n = 3) density 
matrices. It is clearly easily extendible to cases n > 3. The fundamental idea is to scale the 
Oj(f-diagonal elements {pij,i 7^ j) of the density matrix in terms of the square roots of the 
diagonal entries {pu). That is, we set (introducing the new [Bloore] variables Zij)^ 

Pij ^ PiiPjj ^ij ■ ( ) 

This allows the determinant of p (and analogously all its principal minors) to be expressible 
as the product {\p\ = AB) of two factors, one {A) of which is itself simply the product of 
(positive) diagonal entries (pa) and the other — in the n = 4 case under investigation 
here (easily extendible from the case of real density matrices to complex ones) — 

^ ~ (^34 ~ ^12 + 2 (^14 (2^24 — + 2:13 (^23 " 2:242:34)) Zi2 — ^23 — ^24 — 2I4+ (2) 
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^14 (^23 ~ + ^13 (^24 ~ + 22;232;24^34 + 2^132:14 (2:34 — 2;232;24) + 1, 



involving (only) the Zi/s {i ^ j) [20|, eqs. (15), (17)]. Since, clearly, the factor A is positive 
in all nondegenerate cases {pu > 0), one can — by only analyzing B — essentially ignore the 
diagonal entries (and thus reduce by {n — 1)) the dimensionality of the problem of finding 
nonnegativity conditions to impose on p. This is the feature we will seek to maximally 
exploit here. 

It is, of course, necessary and sufficient for p to serve as a density matrix (that is, an 
Hermitian, nonnegative definite, trace one matrix) that all its principal minors be nonneg- 



ative [3l|. The condition — quite natural in the Bloore parameterization — that all the 
principal 2x2 minors be nonnegative requires simply that —l<Zij<l,i^j. The joint 
conditions that all the principal minors be nonnegative are not as readily apparent. But for 
the 9-dimensional real case n = 4 — that is, Q{pij) = — we have been able to obtain one 
such set, using the Mathematica implementation of the cylindrical algorithm decomposition 
3^ . (The set of solutions of any system of real algebraic equations and inequalities can 
be decomposed into a finite number of "cylindrical" parts 33|.) Applying it, we were able 
to express the conditions that an arbitrary 9-dimensional 4x4 real density matrix p must 
fulfill. These took the form, Zi2, zi^, Zi4 G [—1, 1] and 

2^23 £ [^23' ^23]' ^24 G [^24,^24]) ^34 ^ [^34,^3^], (3) 

where 



Z23 = Z12Z1S ±\ 1- zl2\ I- zfr^, Z^^ = zuzu ±\ I- zj2\n- zj^, (4) 



Z: 



± _ Z13Z14 — Zi2ZiiZ23 — Z12Z1SZ24 + Z23Z24: ± S 

Zi2 



and 

S = \J -I + zl^ + 2^3 - 2Z12Z13Z2?, + z\.r^\j-\ + z\^ + z\^ - 2Zi2ZiaZ2A + 2^24- (5) 

Making use of these results, we were able to confirm via exact symbolic integrations, the 
(formally demonstrated) result of Zyczkowski and Sommers 0^ that the HS volume of the 
real two-qubit (n = 4) states is ~ 0.0016106. (This result was also achievable through 
a somewhat different Mathematica computation, using the implicit integration feature first 
introduced in version 5.1. That is, the only integration limits employed were that Zij G 
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[—1, ^ j — broader than those in Q — while the Boolean constraints were imposed 
that the determinant of p and one [all that is needed to ensure nonnegativity] of its principal 
3x3 minors be nonnegative.) 

However, when we tried to combine these integration limits with the (Peres-Horodecki 
n = 4) separability constraint that the determinant (C = IpptI) of the partial 
transpose of p be nonnegative |3Z, Thm. 5] , we exceeded the memory availabilities of our 
workstations. In general, the term C — unlike the earlier term B — unavoidably involves 
the diagonal entries {pa), so the dimension of the accompanying integration problems must 
increase — in the 9-dimensional real n = 4 case from 6 to 9. 



A. Restricting diagonal entries 

Nevertheless, we found that by imposing the condition that the four diagonal entries 
(pii, i = 1, . . . , 4) fall into two equal pairs (say, pu = P22 and ps3 = p^i, so that p^^ = (^-^/^n) ^^ 
the determinant (C) of the partial transpose could now be expressed as the product 

C = \ppt\ = i(l - 2pn)Vn^, < pn < \, (6) 

where 

D = (^24 " 1) ^13 + 22^23 (2^34 — ^14^24) Zis — ^34 Zu~^ (7) 
(Zl4 - 1) (214 + 1) (2:23 - 1) (2:23 + 1) + 2^142:24234 + 

(2I4 — 1) 2:12 + 2 (213^14 + 2:232:24 — (214223 + 2:132:24) 234) 2:12. 

Thus, the term D (like i? in (j2))) is itself independent of the diagonal entries of p (in particular, 
the specific value of pn) — allowing us to proceed, as indicated, with integrations in a lower- 
dimensional {d < 7) setting. (This same form of factorizability takes place, as well, for all 
the [four] 3x3 and [six] 2x2 minors of ppT-) If we can further guarantee the nonnegativity 
of D — in addition to that of i? — we can ensure separability of p. Let us also note that 

B - D = 2(2:14 - 223) (213 - 2:24) (2:12 - 2:34). (8) 

(So, if any of the three factors in (jH)) are zero, the associated state must be separable.) 

We will now proceed to some specific analyses within this more restrictive framework 
(7- dimensional in nature, since we started with the 9-dimensional real setting and have 
essentially only one free diagonal parameter [pn] left). 
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III. 7-DIMENSIONAL REAL SETTING (pn = P22,P33 = pu) 



A. 7-dimensional analysis 

The associated 7-dimensional HS volume of all these states (separable and nonseparable) 

2 

is Y^j2o ~ 0.000652752. The six principal 2x2 minors of the partial transpose simply yield 
the (Bloore) conditions that Zij G [—1, 1], so nothing can be gained — in terms of obtaining 
upper bounds on the separable volume — by using the nonnegativity of these six minors as 
further constraints in our integrations. However, if we require that one of the 3x3 principal 
minors of the partial transpose be nonnegative, we do succeed in obtaining a nontrivial 
upper bound of ^.^2032 ~ 0.000566227 on the 7-dimensional volume of separable states. So, 
we have a derived upper bound (probably rather weak, we surmise) on the HS separability 
probability for our 7-dimensional real set of 4 x 4 density matrices of ~ 0.867446. 

B. 4-dimensional analyses 

1. Z12 = Z23 = Z24 = 

Here, we set the three indicated Bloore parameters to zero. (So, the fourhee parameters 
of the initial seven are pn, 213,^14 and Z34.) Then, using the implicit integration feature 
of Mathematica (rather than the limits (jHl)), we were able to obtain for the total HS vol- 
ume ^ ~ 0.0257021 and (further adding the separability constraint that D, given by ((7j), 
be greater than 0) the separable volume of ~ 0.00902969. Taking the appropriate 
quotient, we find for the HS separability probability for this scenario (quite elegantly), 
^ ^ 0.351321. 

2. Z23 = 0, 2:24 = 0, 2:34 = 

Now, HSl°^i = ifi ~ 0.0218166 and HS'lf = w ~ 0.00902969, so 

HSsepprob — jjoiot ^ 327^- ^ 0.413891. 

vol 

There are twenty possible 4-dimensional scenarios. Of these four each have one of these 
two nontrivial HS separability probabilities (^^ or ^^3^^ '* )• For the remaining twelve cases, 
the HS separability probabilities are simply 1. 
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C. 3-dimensional analyses 



We also observe that for the fifteen possible 3-dimensional scenarios, all the HS separa- 
bility probabilities are (trivially) 1. (Eight of these scenarios have HS total volume equal to 
Y^, four, ^ and three ^.) 

D. 5-dimensional analyses 

For all twelve possible scenarios, setting two Zij^s to zero, the total HS volume of states 
is ifio ~ 0.00685389. 

1- Z23 = 2124 = 

To compute the separable HS volume we had to resort to numerical means, obtaining 
0.00532303, for a separability probability of 0.776643. (This is also the probability for the 
5-dimensional scenario Zi2 = 213 = — and [at least] five others.) 

2. zu = = Z23 = 

2 

Here, as indicated, the total HS volume is ~ 0.00685389. This is also the separable 
volume — since B = D in this case. (This is also the situtation with the scenarios Z13 = 
= and Z12 = z^ = 0.) 

IV. 9-DIMENSIONAL REAL CASE 

As previously mentioned, we know from the Sommers-Zyczkowski analyses P that the HS 
volume of the 9-dimensional convex set of real Ax A density matrices is ~ 0.0016106. (In 
this section — to fully accord with their results 1^ — we have to adjust by an overall scaling 
factor of 2^ = 16 the results of our usual integration procedure employed previously in this 
study, in which we simply employed 1 as our integrand, rather than some other constant. 
This scaling, of course, does not have any impact on separability probabilities.) We can 
also computationally verify the Sommers-Zyczkowski HS volume formula in either of two 
manners: (a) employing the integration limits obtained by application of the cylindrical 
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decomposition algorithm or (b) implementing the implicit integration feature (requiring 
here that the determinant of p and one of its 3 x 3 principal minors be nonnegative) of 
Mathematica, using for the integration limits simply that — 1 < -Zjj < 1 for all i ^ j. We 
investigated how far we could proceed in the full 9-dimensional n = A real case, by imposing 



increasingly greater requirements (corresponding to the Peres-Horodecki criterion |34l l35|) 
that would need to be fulfilled for ppT to be nonnegative definite. Now, of the six 2x2 
principal minors of ppT, only two are distinct from the six such minors of p itself. If we 
demand that, in addition, to the feasibility constraints on p that one of these two minors 
be nonnegative we obtain (to high precision) 0.0014242052589 yielding an upper hound on 
the HS separability probability of the 9-dimensional real two-qubit states of 0.88426997055. 
Further (tighter) numerical results — using additional principal minors of ppx (even the 
remaining 2x2 nonredundant one) — proved difficult to achieve, however. 



V. 3-DIMENSIONAL WEIGHT FOR THE FULL 15-DIMENSIONAL PROBLEM 



The direct straightforward ("brute force") computations by symbolic means of the 15- 
dimensional separable volumes (for any of a wide variety of metrics) of the (complex) two- 
qubit states seem to far exceed present workstation capabilities. Nevertheless, we will here 
try to gain some analytical insight into these formidable problems. 

Let us first note that the computation of the separable volumes could be seen to require 
the evaluation (in which the Peres-Horodecki separability and feasibility criteria are both 
enforced) of a 15-fold integral. Following the innovative Euler-angle-based parameterization 
of 4 X 4 density matrices of Tilma, Byrd and Sudarshan , the first twelve variables over 
which to integrate can be taken to be the twelve Euler angles (oj, z = 1, 12) and the last 
three to be the eigenvalues Ai, A2, A3 of the density matrix p. (Of course, S^^^Aj = 1.) Now 
for any of the metrics of interest, the associated 15- dimensional integrand can be represented 
as the product of a 12-dimensional Haar measure {pHaar) 23|, eq. (34)] {common to all the 
metrics of interest) over the twelve Euler angles 

fJ^Haar = cos(a4)^ cos(a6) cos(aio) sin(2a2) sin(a4) sin(Q;6)^ sin(2a8) sin(a;io)'^ sin(2ai2)(iai2...(iQ;i 

(9) 



(0 < aeven < f,0 < Oodd < 7^ |22 



{pmetric) over the eigenvalues (cf. 



eq. (47)]) and a 3-dimensional metric- specific measure 
38|)- (So, if the Peres-Horodecki criterion is not enforced. 
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we simply obtain the [known in the Bures and HS cases] volumes of the separable and 
nonseparable states, the volumes decomposing into the products of the results of 3-fold and 
12-fold integrations js!, eq. (3.7)], H eq. (3.17)].) 

Consequently, after the first twelve steps of the (presumptively theoretically achieveable, 
but apparently totally impractical) integrations for the separable volumes, we can imag- 
ine obtaining t/iree- dimensional integration problems (over the three-dimensional simplex of 
eigenvalues). The integrands of the problem would now be the products of ^metric and a 
common three-dimensional weighting function W, acquired during the course of the 12-fold 
integration. Certainly, W should be a symmetric function of the four eigenvalues. 



A. First analysis 

We will try to fit W to our various conjectures for the separable volumes [3,0, previously 
obtained by numerical methods. In particular, we have found (after some limited trial 
and error) that the (symmetric) choice (being neither a convex nor a concave function, we 
observed — nor even approximately so). 



53 



W{Xi, A2, A3, A4) = 6086.(AiA2A3 + A1A2A4 + A1A3A4 + A2A3A4)^, (10) 

reproduces the conjectures for both the Hilbert- Schmidt volumes and hyperareas of the sep- 
arable two-qubit states to good accuracy (0.01643%). (These conjectured volume and hy- 
perarea are (5^3)-^ ^ 2.73707- 10"^ and (325^)"^ ^ 7.11111 ■ 10-^ respectively ^, eqs. 
(41), (42)].) For the computation of the 14- dimensional separable hyperarea, the weighting 
function (fTUj) reduces (since we can take A4 = 0) to 

53 

W{X,, A2, A3) = 6086(AiA2A3)5o. (11) 

Now, the acid test of the legitimacy/validity of our choice of weighting function — and the 
raison d'etre of our exercise — is to see how well (in addition to the separable Hilbert-Schmidt 
volumes and h ype rareas (which we constructed to satisfy the Szarek-Bengtsson-Zyczkowski 
two-iold ratio juj]) it reproduces (our presumptively correct) conjectures for metrics other 
than the Hilbert-Schmidt one. 

For the Bures (minimal monotone metric, we found that the use of the weighting 
function (fTO|) — coupled with fisures — predicted a separable volume 0.938275 times the mag- 
nitude of the conjectured value of 2^^^( ^^^ ) ^ 4.2136 ■ 10~^. For the Kubo-Mori monotone 
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metric, we obtained an estimate that is 0.910768 times the magnitude of the conjectured 
value of 2-15(10(72 - 1)) ^ 0.000126408, for the (arithmetic) average monotone metric, 
0.903281 times the magnitude of the conjectured value of 2-^^{f{V2-l)) ^ 0.0000407314, 
and for the Wigner-Yanase monotone metric, 0.919585 times the conjectured value of 
2-^^(1(^/2 - 1)) ^ 0.000221214. So, our choice of W works rather well, at least for a 
first simply heuristic effort. (In [7J, we also had additional volume conjectures for the GKS 
(Grosse-Krattenthaler-Slater) ( "quasi-Bures" ) monotone metric (cf. [41]). However, we en- 
countered numerical difficulties in trying to analyze it here, in the fashion of the other 
metrics.) 

If we try, as well, to predict the 14-dimensional separable hyperarea for the Bures metric 
using (|TT|l . we obtain an estimate of 0.0000262122, which is 0.940364 times as large as the 
conjectured value of2-^^(^^^^§^) ^ 0.0000278746 7, Table VI]. (For the Kubo-Mori metric, 
we also obtain an estimate of the separable hyperarea of 0.0000399861 and of the separable 
probability of a state on the 14-dimensional boundary of 0.0214689 — but there were no 
prior conjectures for these quantities in [3, Table VI]. For the arithmetic average metric, our 
estimate of the separable hyperarea is 0.738784, while our conjecture in Table VI] [to be 
corrected by a factor of 8 as noted in 0, p. 1-11]] amounts to 2-i4( 255(72-1) .^ _ 0.825191.) 



B. Further analyses — exact weighting functions 

We, then, altered our analytical strategy somewhat and succeeded (twice) in reproducing 
our five indicators — - the Bures, Kubo-Mori, Wigner-Yanase and (arithmetic) average sep- 
arable volumes and the Bures hyperarea — all now within 5% of their conjectured values. 
We accomplished this by exactly fitting (by finding the values for a and b) the conjectured 
HS separable volume and separable hyperarea to weighting functions of the form (and their 
A4 = reductions), 

W(Xu A2, A3, A4) = a(Sf<^.A,A,)™^ + 6(Sf<^.<,A,A,Afc)"^^ (12) 

We conducted separate analyses for pairs of low integral values of the exponents (1 < 
"^1, 1^2 < 4). For mi = 3, m2 = 3 we had the results 

a = « 1^26422, * = ^°~ff^ « 19673.7 (13) 

4643757r6 46437571^ ^ ^ 
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^^), and for mi = 4,1712 = 3, 

, 8834477652V3 ^ ^ 33503284082268V3 ^ 

31093757r6 ' 3 1 09 3 757r6 ^ ^ 

(Numerical tests showed that neither of these two functions was Schur-convex 



^ a 3 

nor Schur-concave 4^, nor even approximately so.) As indicated, for both these settings of 
(I12|l . our five (non-HS) indicators all lay within at most 5% of their conjectured values Ql- 
(Additionally, the estimated arithmetic average separable hyperarea lay within 5% of the 
conjectured value for the mi = 4, m2 = 3 case, and 10% for the other. Since the two functions 
both fit the conjectured separable HS volume and hyperarea, any linear combination of them 
will also. We have found that by weighting the function associated with p3|) by 0.570347, 
and the other function by 0.429653, there is no deviation in the associated non-HS indicators 
by more than 3.61%.) The closeness of our estimates to their conjectured values, certainly 
it would seem, should lend some further support (beyond the original numerical evidence 



for the reasonableness of the associated conjectures. 
Of course, it behooves us and possibly other interested researchers, at this stage, to 
explore the properties (and desirably derive guiding principles) of additional candidates 
for the presumptive three-dimensional weighting function W. (We are compelled to note, 
however, that contrary to our construction of W so far, that it is clear that W has to be 



simple flat in some finite neighborhood of the fully mixed state Ai = A2 = A3 = A4 = ;| 
^, sec. 15.5]. So, it seems that W needs to be defined in a piecemeal manner over differing 
domains (cf. 



VI. 5-DIMENSIONAL QUBIT-gt/TO/T WEIGHTING FUNCTION 

Proceeding analogously to our (first) analysis in sec. IV Al for the two-qubit case [n = 
4), we sought to obtain a 5-dimensional weighting function Wn=6 for the qubit-qutrit case 
(n = 6). Using the conjectures — based on extensive numerical results — stated in 0, 
sec. VI. D. 2], in particular, those for the Hilbert- Schmidt separable volume and separable 
hyperarea of the 35-dimensional convex set of 6 x 6 density matrices (conjectured to be 
(2^5 ■ 3 ■ 5^3 . 7^)-i ~ 2.02423 ■ lO'^^ and (2^^ ■ 3 • 5^^y^ ^ 1.94026 ■ 10-^3, respectively), we 
fitted the weighting function (reproducing the HS separable volume to very high accuracy 
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and the hyperarea to an accuracy of .7%), 

W^n=6(Ai, A2, A3, A4, A5, As) = 986304. (S,6<^.<,<;<„A, A,- AfcAA™)i (15) 

Now, in applying this to the Bures metric, we derived an estimate of the 35-dimensional 
Bures volume that was 1.82587 times as large as the conjectured value of 2~^^-3v86429867r ~ 
1.03447- 10~^^ 0, eq. (32)]. Our estimate of the 34-dimensional Bures hyperarea was 1.91223 
times as large as the conjectured value of 2^^ . 3 . 5V84629867r ^ 1.45449 ■ 10"^^ 0, eq. (33)]. 
Though, somewhat disappointingly large, these two early results are certainly of the same 
order of magnitude as the conjectures (and the underlying supporting numerical evidence), 
and suggest additional research. 



VII. REMARKS 



We have sought to determine a certain 3-dimensional weighting function by fitting the con- 
jectured values of the Hilbert-Schmidt separable 15-dimensional volume and 14- dimensional 
hyperarea 0, eqs. (41), (42)]. It would be of interest to attempt to fit additional coYi]ectvjQd. 
values as well (such as those for the Bures [minimal monotone] metric). (No conjecture is 
presently available for the HS separable hyperarea of the 11-dimensional space spanned by 
the rank-2 density matrices. Otherwise it could be incorporated into our further analyses 
too — if the weighting function did not degenerate with two zero eigenvalues present, and 
if the corresponding 11-dimensional separable hyperarea is not actually zero [cf. j^] p, sec. 
VI.C.4]].) 

n 

jet us — as was done in [10|] — bring to the reader's attention some other studies, such 

3, Q, Q, Q, Q 

tainmg to volumes of sets of separable and/or positive-partial- 
transpose states, as well as our more recent analysis [3^, concerning the Hilbert-Schmidt 
metric. (It becomes quite clear in this last study, that the separability constraint on two- 
qubit systems is, in general, quartic in nature (cf. 0,1^]), thus, to some extent, explaining 
the associated difficulties in enforcing it — as wwell as raising certain interesting topological 
questions (cf. 0|)-) 



In conclusion, let us also make reference to a certain capsule review 114 ] in the database 
MathSciNet of our previous paper [7|] in this journal. In particular, we add emphasis to 
the final sentence of the review, devising a response to which comment has been the main 
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motivation of this paper, as well as that of [30|. 

"The paper concerns properties of the convex set of separable two-qubit states. Although 
the positive partial transpose criterion gives in this very case a concrete answer to the 
question of whether a given mixed state is separable, the geometry of the 15-dimensional 
set S of separable states is still not well understood. 

The author analyzes numerically the volume of the set S with respect to measures induced 
by several monotone metrics. In particular, he studies the one-parameter family of metrics 
interpolating between the maximal and the minimal (Bures) metrics. 

Working with the Bures measure he conjectures that the relative volume of the set of 
separable states is equal to the silver mean, a = \/2 — 1. In a similar way the volume of the 
14-dimensional hyperarea of S is estimated with respect to various measures, and the ratios 
area/volume are analyzed. 

The conjectures of Slater, based on numerical integration, still await analytical confirma- 
tion." 
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